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The process of electron-positron annihilation into a pair of tt" and uj mesons is considered in 
the framework of the SU(2)xSU(2) Nambu-Jona-Lasinio model. Contributions of intermediate 
photons, p(770) and p'(1450) vector mesons are taken into account. It is shown that the bulk of the 
cross section at energies below 2 GeV is provided by the process with intermediate p'(1450) state. 
The contribution due to single photon and p(770) exchange is in agreement with the vector meson 
dominance model. Numerical results are compared with experimental data. 
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I. INTRODUCTION 

Studies of the process of associated production of tt" 
and Lu mesons at colliding electron-positron beams pro- 
vide interesting information about meson interactions at 
low energies. Moreover this channel is one of the con- 
tributions to the total cross section of e+e~ annihilation 
into hadrons, which is required for a precise determina- 
tion of the hadronic vacuum polarization. 

The annihilation into the wvr" pair at energies below 
2 GeV was studied experimentally at DM2 ND [3|, 
SND [3], and CMD-2 [4]. The same interactions can be 
also found in the tau lepton decay r — >■ TTuiyr studied at 
CLEO II ll. 

For theoretical description of the process under con- 
sideration the vector-dominance-like models were used, 
see e.g. Ref. To fit the experimental data a set of 
additional parameters describing contributions of ampli- 
tudes with virtual p{770), p'{lA50) and p"(1700) mesons 
was introduced. The energy dependence of these param- 
eters was neglected. Earlier the process of p' — > utt decay 
was considered within a relativistically generalized quark 
model in Ref. 1^6,] and in a non-relativistic quark model [7] . 



In Ref. I8| the reaction 



was considered in 



the vicinity of meson mass region, where the KLOE 
experimental data is available [9]. In this paper we will 
not work specially at this resonance, so that the region 
from the threshold up to about 2 GeV c.m.s. energy will 
be considered without taking into account the </> meson 
contribution. Recently in Ref. |10| the process was con- 
sidered in frames of a non-relativistic quark model. It is 
argued there that the process at energies below 2 GeV 
is dominated by the two-step process in which the pri- 
mary quark-antiquark pair forms a p meson in the ground 
or excited state and then the vector meson decays into 
LU aird TT. It is important to note that the studies in 
papers [1, concluded that the contribution of the 

p"(1700) to the process is small. Following the results 
these works we will neglect the contribution of the am- 




Figure 1: The Feynman diagram with photon exchange. 



plitude with intermediate p"(1700) meson. Meanwhile in 
Ref. it is claimed that for a simultaneous description 
of a series of different annihilation and decay processes 
all three rho meson states should be taken into account. 

In the present paper for the description of the process 
e"'"e~ — >■ wtt" we will use the version of the Nambu-Jona- 
Lasinio (NJL) model, which allows us to describe both 
the ground and the first radial-excited mesoir states 
U^. Note that for the description of the amplitudes 
with virtual photon and the ground p(770) state one can 
use the standard NJL model [171424] . It worth to note 
that for the case of the ground meson states both ver- 
sions of the NJL model lead to the same results, see e.g. 
Refs. In our model it is possible to describe as 

the transition amplitudes 7* — > p, p' as well as the ver- 
texes j*,p,p' — > TT^w without introduction of any addi- 
tional arbitrary parameters. Moreover, the description of 
the vertexes using quark triangle diagram of the anomaly 
type allows us to get the energy dependence of them. 



II. PROCESS AMPLITUDES 

For the description of the first two diagrams, see Figs.[T] 
and [21 we need the part of the standard NJL Lagrangian 
which describes interactions of photons, pions and vector 
p and uj mesons with quarks, see Refs. (itI. [21|. [2^ . It has 
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Figure 2: Feynman diagram with p meson exchange. 




Figure 3: Feynman diagram with p meson exchange. 



the form 
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eQA + igT,j5T3Tr° 



T3P 



(1) 



where q = {u,d) with u and d quark fields; m = 
diag{mmrnd), rriu = rrid ~ 280 MeV is the constituent 
quark mass; Q = diag{2/3, —1/3) is the electromagnetic 
quark charge matrix; e is the electron charge; A, tt*^, lo 
and p*' are the photon, pion, w and p meson fields, re- 
spectively; (77r is the pion coupling constant, = niu/ f-K, 
where = 93 MeV is the pion decay constant; gp is the 
vector meson coupling constant, gp « 6.14 corresponding 
to the standard relation (7^/(47r) = 3; / = diag{l, 1) and 
Ts is the third Pauli matrix. 

All three amplitudes contain the common part corre- 
sponding to e^e^7 vertex and the photon propagator. 
So the sum of the amplitudes can be cast in the form 



(2) 



where s — (pi(e+) -f p2(e^))^ = q^- The first part Ti 
is just the triangle quark diagram of the anomaly type. 
Note that the loop integral in it is finite. Following pa- 
pers [25-27,], we use here the naive confinement approach 
and neglect the imaginary part of the loop integral. 

The integral over the energy fc" of the virtual loop 
momentum is calculated analytically using the residue 
method. The integral over k is taken numerically. Even 
so that this integral is convergent, we put the cut off for 
the upper value of |^| equal A = 1.03 GeV [H. This 
cut off will be necessary in contributions of the radially 
excited mesons states. Here the cut off is applied for 



homogeneity of the approach. And the numerical result 
for the convergent integral do not change considerably 
if the cut off would be removed. The imaginary part is 
neglected by taking the principal value of the integral. 

The second contribution T2 contains three factors. The 
first one is the transition of photon into p meson which 
is described in Ref. (T7| : 
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Note that contrary to the case of the triangle diagram, 
the quark loop describing the 7 — p transition contains a 
logarithmic divergence. The standard NJL methods were 
applied for its regularization using the cut off value. The 
second factor is the p meson propagator. 



ig- 



M2 + iMpTp 



(4) 



where the neutral p-meson mass Mp = 775 MeV and 
width Tp = 146 MeV [Hi. Note that the non-diagonal 
terms in the numerator of the vector particle propagator 
were dropped because of the gradient invariance of the 
triangle diagram. The third factor is the same triangle 
diagram as in the first amplitude Ti. 

A more complicated situation appear for the third con- 
tributions T3, see Fig. [31 because we deal here with ra- 
dially excited p' meson. Instead of the Lagrangian ([T]) 
we use here an extended version of the NJL Lagrangian 
which allows us to describe both ground and radial- 
excited meson states [13, [l3, \2M: 
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q{k')\ id-m + eQA + A^T^j5Tr°{p)+A^Lo{p) 



(5) 



sin(a-l-Q!o) , ^2 sin(a-Q;o) 
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The radially-excited states were introduced in the NJL 
model with the help of the form factor in the quark-meson 
interaction: 

fik^') = {l-d\k^'\)Q{A'-\k^'\), 



k^ = k- 



{kp)p 



d = 1.78GeV~^ (6) 



where k and p are the quark and meson momenta, respec- 
tively. The filled circles in Fig. [3] denote the presence of 
the form factor in the quark- meson vertexes. Note that 
the NJL model itself and its extended version can be 
used only for sufficiently low energies. In this study we 
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attempt to receive qualitative results working at energies 
up to 2 GeV. 

Coupling constants 17^^ and gp-^ coincide with g^^ and gp 
constants introduced above in the standard NJL version. 
The other coupling constants are defined via one-loop 
integrals: 



9^2 = 

where 



4/, 



-1/2 
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-1-1/2 
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d^fc (/(fc-^'))" 
(27r)4 (to2 - /fc2)m ' 



= \/65^2, (7) 



•n,m — \, 2. 



The angles ao = 59.06°, a = 59.38°, /3o = 61.53° and /3 = 
76.78° were defined in Ref. [3 [13 to describe mixing of 
the ground and excited meson states. This contribution 
Ta again consists of 3 parts. The 7-/32 transition (7 — pi 
transition coincides with the standard 7 — p one) can be 
expressed via the 7 — p transition ([3]) with the additional 
factor [Hill 



r = 



0.47. 



(8) 



So the 7 — transition takes the form 



— {9 q ~q q 
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We take the p' propagator is taken in the Breit-Wigner 
form 



9^ 



(9) 



where the running p' width reads 

rp,{q^) = r{p' ^ 27t) + r{p' ^ + {Tp,{Mp,) 



PaAMp')J' 

where pa, (s) is the momentum of oi meson in the decay 
p' —T' aiTT. We assume that below the threshold of the 
reaction p' — >■ ain the main contribution the the width is 
given by the two channels p' — ?> 27r p' — >■ wtt". Above the 
peak ^/s > Mpi , where many other channels are opened, 
we use the complete width T pi = 340 MeV (we take the 
value at the lower PDG [2^ boundary). The transition 
to the complete width is approximately described by lin- 
ear switching on of the contribution due to the decay 
p' — aiTT being one of the most probable channels. The 
values r(p' ^ 27r) = 22 MeV and T{p' ^ wtt^) = 75 MeV 
were calculated in |14] in agreement with the experimen- 
tal data [S^l ■ Since we are working close to the ujit thresh- 
old, taking it into account in the running width is impor- 
tant. Running of the p meson width is less important 
numerically, since the p meson contribution is relatively 
small. 



III. NUMERICAL RESULTS AND DISCUSSION 

Now we can estimate the contributions of the consid- 
ered amplitudes into the total process cross section. The 
details of phase volume calculations and evaluation of the 
cross section can be found in Ref. [l^j- For our case it 
takes the form 



xBr(w 7r°7), 
\{s,MlMl)^{s 



(11) 
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where 



J(3) = ( 1 _ 



g2 - M 2 + iMpTp I ^ 
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q'-^e+^^Tp,iq^) ' 
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and 



/(3) 



d'^fcm2e(A^ - \k^ \] 
i7r2 (fc2 — m2 + iO) 
1 



{{k - P^f - m2 + iQ){{q - kY -m? + iO) 

In the first line of Eq. <\V2\ we have the sum of the pho- 
ton and rho meson exchange contributions. Their sum 
takes the form that coincides with the one received in 
the vector meson dominance model, see e.g. In fact, 
the standard NJL model contains the vector dominance 
model [17,,|3l|,|33. 

Note that keeping the cut-off for the convergent inte- 
gral in I\ entering Ti and T2 in Eq. ^ is not neces- 
sary, but it does not affect much the numerical result. 
Expression for the integral has a rather cumbersome 
form and contains a combination of terms with different 
powers of the form factor (up to the third power). It is 
constructed according to the Feynman rules coming from 
the Lagrangian ([5]) . For calculation of the relevant quark 
loop integrals we use the method described in Ref. (25| . 
It is worth to note that in our calculations the signs of 



I'f^ and appeared to be opposite in accordance with 
the fit to experimental data performed in [4]. 

The coupling constants gp — Q and /^r — 93 MeV 
entering in Eq. [TT] are universal input parameters for 
the NJL model. In Ref. [4] another value for this con- 
stant was used: /p w 5 received from the decay width 
r(p — > e^e~). Another difference is coming from the 
value for the coupling constant in the vertex pcjTr. In our 
model it is gp^^, = 35ip/(87r2/^) « 14.7 GeV"\ while in 
Ref. ^] the value gp^T^ « 17 GeV-i taken as a fitting 
parameter. 
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Figure 4: Comparison of experimental results for e^e 
tt'^oj — >■ 7r''7r"7 with the NJL model prediction (lines). 



Fig. m shows the experimental data P, H, 01 and the 
corresponding theoretical prediction (the solid line) re- 
ceive within the applied here NJL phenoinenological 
model. The dash-dotted hire shows the sum of the photon 
and rho-meson exchange contributions. The short-dash- 
dotted corresponds to the pure p' meson exchange. The 
photon and p meson exchange is important for the thresh- 
old region, while the p' contribution dominates in the re- 
gion \/s ^ Mpi . Note that the NJL model is adjusted for 
applications at low energies up to about 2 GeV. In this 



energy range, the model gives a qualitative description 
of meson properties and interactions. The advantage is 
that the set of parameters is limited and fixed. Note that 
to describe the given process we did not introduce any 
new parameter in the model. Presumably, adding of the 
p"(1700) meson contribution might improve the agree- 
ment with the experimental data above the peak, but for 
the time being the NJL model is not suited to include the 
second radial excitations of mesons with large masses. 
A more accurate description of the threshold behavior 
requires going beyond the Hartree-Fock approximation 
that was used here. Indeed, meson-meson final state in- 
teractions can play an important role in the threshold 
domain. 

The same approach was successfully applied in pa- 
pers for description of mass spectra and strong de- 
cays with participation of excited mesons. In the present 
work we continue the work started in Refs. (25l . [33| de- 
voted to description of radiative decays with participa- 
tion of radially excited mesons and pass to description of 
annihilation processes studied at modern e^e~ colliders. 
Similar mechanism appear in the processes of e^e~ an- 
nihilation into e.g. 7r°7, 7r'7, and tt'^tt^ which will be 
considered elsewhere. 
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